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In this work we present a theoretical framework that aUows for the existence of coherent twinlike 
models in the context of self-dual Maxwell-Higgs theories. We verify the consistence of this frame- 
work by using it to develop some twinlike self-dual Maxwell- Ifiggs models. We use a combination of 
theoretical and numerical techniques to show that these models exhibit the very same topological 
BPS structures, including their field configurations and total energy. The study shows that it is 
possible to develop a completely consistent prescription, which extends the idea of twinlike models 
to the case of vortices in Maxwell-Higgs theories. 

PACS numbers: ll.lO.Kk, ll.lO.Lm 



I. INTRODUCTION 

Topological structures have been extensively used to 
study phenomena related to several areas of physics [l| . 
In particular, they are of direct interest to cosmology, 
since they can be formed in a rather natural way dur- 
ing phase transitions in the early universe. In the con- 
text of classical field theories, topological structures are 
described as finite energy solutions to some nonlinear 
field models. Usually, the models are endowed with a 
scalar potential which exhibits degenerate vacua, allow- 
ing for the spontaneous symmetry breaking mechanism. 
In this context, kinks which are the simplest topo- 
logically non-trivial configurations, are the 1-dimensional 
solutions to a single real scalar field theory, while vortices 
[sj and monopoles p| are, respectively, 2-dimensional so- 
lutions to some Abelian gauge theories and 3-dimensional 
solutions to non- Abelian gauge theories. 

An interesting issue concerning the study of topological 
structures is that, in some particular cases, such struc- 
tures, beyond solving the standard equations of motion, 
are also solutions to a set of coupled first-order differen- 
tial equations, named Bogomol'nyi-Prasad-Sommerfield 
(BPS) equations In this context, since we do not need 
to deal with the second-order Euler-Lagrange equations 
of the model under investigation, the topologically non- 
trivial BPS configurations are easier to be obtained. An- 
other advantage of studying BPS configurations is that 
they describe minimum energy solutions to the corre- 
sponding classical field theories. 

Motivated by recent results on cosmology and super- 
string theories, a new type of classical field theories has 
been intensively investigated during the last years. These 
theories, named k-field models, are usually endowed with 
non-standard kinetic terms that change the dynamics of 
the model under investigation. In this context, topolog- 
ically non-trivial BPS solutions may exist, and the fea- 
tures they engender can be quite different, or quite simi- 
lar to the ones presented by their usual counterparts; see, 
for instance, Ref. [6[. 

In a recent work, Andrews, Lewandowski, Trodden and 



Wesley noted that some k-theories can mimic the very 
same topological structures engendered by their stan- 
dard counterparts, including their energy densities Q. 
Since these models map each other solutions, they were 
named twinlike models. In another recent work, some of 
us studied the existence of twinlike models in the con- 
text of scalar field theories 8]. There, in particular, one 
used a combination of theoretical and numerical tech- 
niques in order to perform a detailed analysis of such 
models, including the presence of first-order differential 
equations and their applications to the braneworld sce- 
nario. Soon after, another work by Adam and Queiruga 
made very good use of the first-order framework, to 
develop algebraic technique to construct twinlike theo- 
ries. Yet more recently, another interesting result was 
introduced in Ref. concerning the presence of twin- 
like models, where the models present not only the same 
solution with the same energy density, but also the very 
same stability features. 



The above works on twinlike models deal with scalar 
fields. Thus, a rather natural question is how to con- 
struct consistent twinlike models under the action of 
gauge fields. In the present work we focus our attention 
on planar Maxwell-Higgs models and their topologically 
non-trivial self-dual solutions. To achieve this goal, we 
organize the work as follows: in the next section we in- 
troduce the generalized theory and we also present the 
theoretical framework which allows for the existence of 
twinlike self-dual Maxwell-Higgs models. In Sec. IIIIl we 
prove the consistency of such framework by using it to de- 
veloped a few twinlike self-dual Maxwell-Higgs systems. 
Next, in Sec. IIVI we perform a detailed numerical analy- 
sis of these systems and in Sec.|Vl we present conclusions 
and some perspectives concerning future investigations. 
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II. THE THEORETICAL FRAMEWORK 

Let us first present the model, which is defined by the 
planar Lagrange density 



hm 



wm)\D,^\'-vM) . (1) 



Here, F^i, — d^j^Ai, — d^A/^i stands for the electromagnetic 
field strength tensor, D^(f> = d^(j) + ieA^(j) represents the 
usual covariant derivative and ^ (|0|) is the spontaneous 
symmetry breaking potential which controls the scalar 
matter self-interaction. Also, ft.(|(/>|) and w are di- 
mensionless functions of the amplitude of the scalar field. 
We note that non-trivial choices to h{\(j)\) and ui (|(/)|) will 
induce, respectively, non-standard dynamics to both the 
gauge and scalar fields. Also, there are interesting mo- 
tivations concerning the use of these two functions and, 
in particular, the function h (|0|) represents a generalized 
dielectric function, which can be used to describe inter- 
actions between quarks and gluons [T]| . 

Since it is useful to deal with dimensionless fields, co- 
ordinates and parameters, let us introduce the mass scale 
M, which we use to implement the following scale trans- 
formations: a;^ -> x^/A/, Af^0, M^A^^, 
e — >■ M^e and v — ^ M^v, where v represents the sponta- 
neous symmetry breaking parameter of the model under 
investigation. As a consequence, we get that £ -^M^Cg, 
where Cg stands for the dimensionless Lagrange den- 
sity to be used from now on, which has the very same 
functional form of C. Also, for simplicity, we choose 
V = e = 1. 

The standard procedure to search for the rotationally 
symmetric configurations is to consider such configura- 
tions as static solutions of the Euler-Lagrange equations 
of the model. In the present case, these equations are 
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where = F^.^Fi"'. Also, we take = -2w |0|^ A^. 

In the present case, the Gauss law for static fields can 
be written in the form (k represents spatial indices) 



(4) 



from which we note that such law is trivially satisfied by 
A^ — (temporal gauge). So, we fix this gauge and use 
it from now on. 

We now look for static and rotationally symmetric con- 
figurations of the form 



ct>{r,e)=g{r) e^' 



(5) 



(a (r) - n) 

r 



(6) 



where r and are the polar coordinates. Also, n = ±1, 
±2, ±3, represents the vorticity of the configuration. 
We use the above eqs. ([5]) and © into the Euler-Lagrange 
equations ^ and ^ to get 
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These are the rotationally symmetric equations of motion 
for the fields a (r) and g (r) , respectively. 

In order to solve the Euler-Lagrange equations ([7]) and 
([8]), we need to specify the model. In general, it can be 
done following the generalized model, i.e., choosing non- 
trivial functional forms to the functions h (g) and w {g). 
In this case, it is important to note that both functions 
must be positive-definite, in order to avoid problems with 
the positiveness of the total energy of the model; see 
the expression for the energy density below. We have 
to choose h{g) and w{g), and we also need to specify 
the potential V (g) for the scalar matter self-interaction, 
which must allow for the spontaneous symmetry breaking 
mechanism. 

Before specifying the model, let us recall that the 
main purpose of this work is to investigate the existence 
of coherent twinlike models in the context of classical 
Maxwell-Higgs field theories. In particular, we focus our 
attention on the study of twinlike BPS states, i.e., fi- 
nite energy field configurations which can be described 
as static solutions of a set of first-order differential equa- 
tions (the BPS equations) having the very same field pro- 
file and energy. 

In this sense, it is instructive to note that the limit 
w (g) = 1 leads the model ([T]) back to the model studied 
in Ref. Thus, the limit h{g) — 1 leads us back to 

the standard Maxwell-Higgs electrodynamics. Thus, a 
very interesting choice to the symmetry breaking poten- 
tial V (g) for the scalar matter self-interaction is 
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Now, according to the conventions stated above, the 
Euler-Lagrange equations of motion ([7]) and (|8]) for the 
rotationally symmetric solutions can be rewritten in the 
form 
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The equations (ITO|) and pT|) are solvable by solutions of 
the two first-order differential equations 



dr 



± 
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boundary conditions to be imposed on the functions g (r) 
and a (r) , near the origin and asymptotically. Near the 
origin, such functions must have no singularity. So, given 
the Ansatz ([5]) and ([H]), g{r) and a (r) have to behave ac- 
cording to 



1 da , / 2 

~:r = =^ (-9 

r dr 



(13) 



The static solutions g (r) and a (r) to the above 
first-order equations are the well-known Bogomol'nyi- 
Prasad-Sommerfield (BPS) states related to the standard 
Maxwell-Higgs electrodynamics. These states solve the 
rotationally symmetric Euler-Lagrange equations of mo- 
tion ([To]) and pT|) by minimizing the total energy of the 
resulting solutions; see equations ([22l) and ([23]) below. 

If we use non-trivial choices for the functions h (g) and 
w (g), the resulting Euler-Lagrange equations ([7]) and ^ 
will be much more sophisticated than the usual Eqs. ([TU]) 
and (ITU) . In this sense, it seems useful to calculate the 
general expression for the non-standard energy density, 
from which we can get interesting insights related to the 
non-usual BPS states. This expression is 
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Here we note that the presence of the functions h (g) and 
w (g) makes it very hard to obtain the BPS configurations 
related to the non-standard model (H]). However, the ex- 
istence of coherent twinlike BPS configurations, i.e., finite 
energy static solutions of the first-order eqs. (fT2l) and ([T3l) . 
is still possible, and it is closely constrained by two im- 
portant assumptions concerning the functions h{g), w (g) 
and V {g). The first one is the following differential rela- 
tion 



= h 
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which leads to a link between the functions h (g) and 
w {g). The second constraint is 



Vig)^h{g)V,{g) 
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which defines the potential V (g) of the non-standard 
model (ITJ as a product between the function h (g) and 
the self-dual potential (O of the canonical Maxwell-Higgs 
model. As we demonstrate below, via the constraints 
(ITS!) and (fTBj) . the non-standard energy functional (fT^ 
can be rewritten in a way such that the corresponding 
total energy is minimized by the first order equations 
(1121) and (|13p . Furthermore, such minimum energy can 
be adjusted to mimic the very same value achieved by the 
standard BPS solutions. This adjustment can be done by 
tuning the boundary conditions to be imposed on h {g), 
near the origin and asymptotically. We perform this tun- 
ing in the next Section. 

Now, to search for twinlike BPS solutions to the first- 
order equations (|12p and (fT5)) . we need to know the 



g (r -> 0) and a (r -> 0) -J> n . 



(17) 



Also, the twinlike BPS configurations must have finite to- 
tal energy. In this case, as a condition to make the energy 
finite, the non-standard energy functional (fT^ must van- 
ish for r — > oo. So, asymptotically, the profile functions 
g (r) and a (r) must behave according to 



g {r ^ oo) — > 1 and a (r — > oo) 







(18) 



III. TWINLIKE BPS MODELS 



Let us now focus our attention on the twinlike BPS 
configurations themselves. The main purpose here is to 
use the theoretical framework developed in the previous 
Section to assure the existence of such twinlike configura- 
tions. We start using (|15l) and ([TS]) to rewrite the energy 
density (fTl)) in a way such that the resulting total energy 
is minimized by the differential eqs. ([T^ and (fT5)) . The 
value achieved by such minimized energy depends on the 
behaviour of h (g), near the origin and asymptotically. In 
this sense, by choosing appropriate boundary conditions, 
we tune this value to be the very same one gets in the 
context of the standard self-dual Maxwell-Higgs model 

(ED. 

In general, since the functions h {g) and w (g) are re- 
lated to each other by the differential constraint ([T5|) , and 
the non-standard potential V (g) is given by Eq. (|16p . the 
modified energy functional (jl4l) can be rewritten in the 
form (see also Eq. (jH])) 



h f 1 da / , 
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do qa 
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from which we note that the corresponding total energy 
E is minimized by the first order differential equations 
(1121) and (I13p . In this case, the minimum energy is given 
by 



bps 



where 



(20) 



(21) 



is the minimum energy density, i.e., the energy density 
related to the first order equations (IT^ and (fT5)) . Also, 
Eq. (1201) defines the Bogomol'nyi bound, i.e., the lower 
bound for the energy functional (jl4[) . 
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We use Eq. (|211) to see that the Bogomornyi bound 
(PO)) depends on the boundary conditions ([T7]) and ([T^. 
and also on the boundary conditions to be imposed on 
h (g) itself. Here, the interesting point is that there 
are infinitely many functional forms to h {g) which obey 
the same general boundary conditions; as a consequence, 
there are infinitely many non-standard energy densities 
(all of them given by Eq. (HI])) which achieve the very 
same Bogomol'nyi bound (given by Eq. ((20)) ). Each one 
of these forms of h (g) defines a different non-standard 
self-dual model. In this sense, since the first order equa- 
tions (IT^ and (fT5|) . and also the boundary conditions 
([TT)) and (fT8|) . do not depend on h (g), all these models 
exhibit the very same BPS solutions to the profile func- 
tions g (r) and a{r). Since these different models have 
the very same field configuration and total energy, we 
keep naming them twinlike models. 

To better clarify how the theoretical framework devel- 
oped in this paper works, let us first consider the stan- 
dard Maxwell- Higgs model, which is defined hy h{g) — 1. 
In such context, due to the finite-energy boundary con- 
ditions dlTl) and dUl), Eq. ([20]) gives the Bogomol'nyi 
bound (which is quantized according to the vorticity n of 
the solutions) 

Ebp, = 27r|n| , (22) 

where 

1 da Id,,, , , 

Ebps = T-:r±-x 5 " ' 23 

r ar r ar 

in the minimum energy functional. Here, we point out 
the existence of infinitely many models (i.e., infinitely 
many functional forms to h{g)) which achieve the very 
same Bogomol'nyi bound (|20)) . In fact, as the reader can 
easily check, any function h (g) satisfying the boundary 
conditions (C is any non-negative constant) 

h{r -^0) and h{r ^ oo) (24) 

leads to E^ps = 27r 

In order to illustrate the above result, let us now con- 
sider some specific models. The first example, represent- 
ing the first non-standard self-dual model, is defined by 

h{g)^g' + l. (25) 

In this case, C = 2. The energy functional (PT]) can be 
rewritten as 

1 da 1 d , A . , 
r ar r ar 

Then, via eqs. (1^01) and I^E^, we conclude that the modi- 
fied model defined by ((25|) achieves the same Bogomol'nyi 
bound (P^ achieved by the usual Maxwell-Higgs model. 

The other non-standard self-dual model we present is 
defined by 

h{g)^{g' + l){g'~lf , (27) 



with C = in this case. Here, the minimum energy 
density is 

e,,. = T;^±i|:(a5M/-2/ + 2)) . (28) 

It is important to note that, due to the boundary con- 
ditions dlT]), dH]) and ([24]), the second term in eqs. ([23]) . 
(1^51) and (1^5)) does not contribute to the integration pro- 
cess in (|20)) . In this sense, we conclude that all the phys- 
ical information concerning the Bogomol'nyi bound (j22p 
achieved by the models previously presented is enclosed 
by the first term in eqs. (^5]) . (P5|) and (P5|) . So, from now 
on we refer to such first and second terms as the physical 
and the irrelevant ones, respectively. In fact, all models 
defined by a /i (g) of the form 

h{g) = l + Gig) (29) 

exhibit a minimum energy functional given by 

Ida Id/,/ q'^ -I W 
^^.^=T-^±-^(«.9^(l + Vjj ■ (30) 

Then, as a consequence, such models possess the very 
same physical term, while exhibiting quite different irrel- 
evant terms. 

To end this Section, we point out that one can choose 
the boundary conditions (j24[) in a different way. Generi- 
cally, one can define {N is any positive constant) 

h{r^O)^N and h{r^oo)^C. (31) 

Here also, there are infinitely many functional forms to 
h (g) which behave according to such boundary condi- 
tions. Each one of these forms defines a different non- 
standard self-dual model. These models compose an in- 
finity set of twinlike theories, since they exhibit the very 
same BPS configurations (given as solutions of the differ- 
ential eqs. (IT^ and (IT51) . according to the finite-energy 
boundary conditions ([TT]) and (fT8]) ) and also achieve the 
very same Bogomol'nyi bound (given by Eq. (pHl)). which 
is 

Ebps^^nlnlN . (32) 

If we use a general N, none of the above models encloses 
the standard Maxwell-Higgs model, which is defined with 
iV = 1. We also note that the extended versions of (|29)) 
and (|30]) are quite obvious. 

IV. TWINLIKE NUMERICAL SOLUTIONS 

Let us now perform the numerical study of the twinlike 
self-dual models presented in this work. We numerically 
solve the first-order differential equations (fT2)) and (fT3]) . 
according to the finite energy boundary conditions (|17p 
and (fT8|). The resulting solutions are well-known, and 
can be found, for instance, in Ref. We then use such 
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solutions to perform a detailed analysis of the energy of 
the vortices, including their minimum energy densities 
([23| . (|26l) and (|28|) . In this sense, we also depict, sepa- 
rately, the numerical solutions to the physical and to the 
irrelevant terms presented in the previous Section, from 
which we show (numerically) that the irrelevant terms 
do not contribute to the Bogomol'nyi bounds computed 
above. Some results are plotted below, for the case n = 1. 

In Fig. 1, we present the numerical solutions for the 
energy densities related to the twinlike BPS states, and 
we note that the solutions for ([^ and engender the 
same basic features: they reach their maximum values 
near the origin, and they decrease monotonically for all r. 
Here, it is interesting to note that the core of the twinlike 
solution (pS)) is smaller than that of the usual solution 
On the other hand, we see that the solution for 
the minimum energy density (I26[) is quite different from 
the previous ones. In this case, such solution reaches its 
maximum value at some finite distance R from the origin, 
and it is not monotonically decrescent for all values of 
the independent variable r. Finally, we point out that all 
the minimum energy densities previously depicted vanish 
asymptotically. 
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FIG. 1: Plots of solutions to the minimum energy densities 
(f23ll (black solid line), ^ (red dashed line) and JMl (blue 
dotted line). 

In fact, since the minimum non-standard energy func- 
tional (j21l) is r-dependent only, the information regarding 
the Bogomol'nyi bound ([20|) is enclosed by the product 
rsbps- In particular, the total energy Ebps of the twin- 
like BPS states is given by 2tt times the area enclosed 
by TEbps ■ Regarding the non-standard models previously 
presented, such area is always equal to |n|. Thus, all 
the above models achieve the very same minimum total 
energy, which is given by the Bogomol'nyi bound (1221) . 



1.2 




r 



FIG. 2: Plots of solutions to the product retps- Conventions 
as in FIG. 1. See also the solution to —a' (green dot-dashed 
line). 



In Fig. 2, we depict the numerical solutions for the 
product rebps related to the previous models. In general, 
all the solutions engender the same features: they reach 
their maximum values at some finite distance from the 
origin, while vanishing for r — !■ and for r — > oo. It is 
interesting to note how these solutions behave in order 
to enclose the same area ever: there is an inverse relation 
between their maximum values and their characteristic 
lengths, and, as a consequence, the solutions with greater 
amplitudes exhibit smaller cores, and vice- versa. 

In fact, the product rsbps means a summation over 
two different terms: the first one is given by r times the 
physical term, and the second one is given by r times 
the irrelevant term; see eas. lj^ . and (1^ . So, in 
order to get to some useful insights concerning the total 
energies of the twinlike BPS states, let us evaluate the 
numerical solutions for such terms, separately. 

In the same Fig. 2, we plot the solution to —a' (prime 
means a derivative with respect to r). This solution is 
important since it denotes the product between r and 
the physical term, and we note that it exhibits the same 
basic features the solutions to rSbps do, including the in- 
verse relation between the amplitude and the characteris- 
tic length. Here, it is interesting to note that this solution 
also encloses an area equal to |n|. So, as stated theoreti- 
cally before, our numerical analysis reinforces that all the 
physical information regarding the total energies of the 
non-standard models previously presented is enclosed by 
the physical term itself. 

In Fig. 3, we depict the solutions to the second term, 
that is, the product between r and the irrelevant term in 
the energy density. In this case, it is interesting to note 
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re. 




FIG. 3: Plots of solutions to the product between r and the 
irrelevant terms. Conventions as in FIG. 1. 



how these solutions behave in order to contribute noth- 
ing to the Bogomol'nyi bound. Here, we reinforce that, 
even in the presence of the negative sectors presented by 
such solutions (and, as a consequence, by the irrelevant 
terms themselves, since the radial coordinate r is always 
positive), the positiveness of the minimum energy func- 
tionals (^5]) . (pSj) and is completely assured; see, for 
instance, the solutions depicted in Fig. 1. 

To end this Section, let us investigate a very important 
issue regarding the study of topological configurations, 
which is the existence of conserved topological charges. 
In fact, the existence of such charges assures the topo- 
logical stability of the corresponding configurations. In 
the present context, that is, for rotationally symmetric 
solutions of the form ^ and ^ , the corresponding topo- 
logical charge Qt can be identified with the flux of the 
magnetic field generated by such solutions themselves. 
So, one has 



Qt = / B(fr = $i 



where $b stands for the magnetic flux, and 

Ida 

B{r) = — — 
r ar 



(33) 



(34) 



is the magnetic field; see ([5]). According to our conven- 
tions, the magnetic flux (and, as a consequence, the topo- 
logical charge) can be written as 



2iTn , 



(35) 



which shows that Qt is conserved, and is quantized ac- 
cording to the vorticity n. Since this result does not 



depend on h (g), we think that the topological stability 
of the non-standard solutions is achieved in the same way 
as that of the standard case. 



V. ENDING COMMENTS 

In this work we investigated the existence of consis- 
tent twinlike self-dual classical field theories in the con- 
text of the planar Maxwell-Higgs model. We have used a 
modified Maxwell-Higgs model endowed by non-standard 
dynamics to both the gauge and the scalar fields. The 
unusual dynamics were introduced in terms of two non- 
trivial functions, /i(|(/)|) and wd^j), which are functions 
of the amplitude of the scalar field; they must be pos- 
itive, in order to assure the positiveness of the energy 
functional of the modified model. We note that ft. (|0|) is 
a kind of generalized dielectric function which couples to 
the canonical Maxwell action, leading to a non-standard 
dynamics to the gauge field. Also, w {\<j)\) couples to the 
squared covariant derivative of the scalar field. 

We have demonstrated the existence of twinlike self- 
dual Abelian-Higgs models under the presence of two 
specific constraints: the first one is a differential rela- 
tion between the non-trivial functions h (|0|) and w {\4>\), 
and the second one is a relation between ft (|0|) and the 
spontaneous symmetry breaking potential (|0|) of the 
non-standard model. We have used these constraints to 
performed a detailed numerical investigation, studying 
the energy spectra of the generalized models, including 
their minimum energy functionals. In general, we have 
noted that the numerical profiles for Sbps itself can be 
quite different from one another, as they can engender 
different basic features. Even in this case, we have no- 
ticed that the solutions to rsbps engender the same fea- 
tures, including an interesting inverse relation between 
their maximum values and their characteristic lengths. 
In particular, such relation is very important, since it as- 
sures that all the solutions to rsbps enclose the same area, 
so all the non-standard BPS states have the very same 
total energy. 

We hope that the above results may stimulate subse- 
quent analysis in the field, mainly regarding the features 
that the twinlike Maxwell-Higgs models may engender. 
In particular, a rather natural issue concerns the exten- 
sion of the present idea to the case of non-rotationally 
symmetric solutions. Another issue concerns the exten- 
sion of the twinlike models to the context of Yang-Mills- 
Higgs theories. These and other related issues are under 
investigation, and we hope to report on them in the near 
future. 
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